Capillary filling in patterned channels 
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We show how the capillary filling of microchannels is aflFected by posts or ridges on the sides 
of the channels. Ridges perpendicular to the flow direction introduce contact line pinning which 
slows, or sometimes prevents, filling; whereas ridges parallel to the fiow provide extra surface which 
may enhances filling. Patterning the microchannel surface with square posts has little effect on the 
ability of a channel to fill for equilibrium contact angle 9e < 30°. For 9^ > 60°, however, even a 
small number of posts can pin the advancing liquid front. 
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Recent years have seen rapid progress in the technol- 
ogy of fabricating channels at micron length scales. Such 
microfluidic systems are increasingly finding applications, 
for example, in diagnostic testing and DNA manipulation 
and as microreactors. As narrower channels are used, 
to conserve space and reagents, surface effects will have 
an increasing influence on the fluid flow. In particular, 
it may be possible to exploit the chemical or geomet- 
rical patterning of a surface to control the flow within 
the channels. Indeed first steps in this direction have 
already been taken: for example, Joseph et al. used su- 
perhydrophobic surfaces to introduce slip lengths of the 
order a few /im [Tj, Zhao et al. used hydrophilic stripes 
to confine fluid flows , and Stroock et al. used oriented 
grooves to enhance mixing Q . 

In this paper we describe the way in which a fluid pen- 
etrates a microchannel patterned with posts or ridges. 
We find that the hysteretic behaviour that results from 
the surface patterning can be a substantial factor in con- 
trolling how the fluid moves down the channel. Ridges 
or posts oriented parallel to the flow speed up capillary 
filling while those perpendicular to the flow pin the in- 
terface and suppress, or sometimes prevent, filling. We 
investigate the rate at which the microchannel fills, find- 
ing that it is highly dependent on both the geometry and 
the spacing of the posts. In particular, we demonstrate 
that a patterned surface can act as a valve, allowing fill- 
ing in one direction only. 

The classical analysis of capillary penetration is due 
to Washburn Consider a capillary of width h with 
an infinite reservoir of liquid of dynamic viscosity 77 at 
one end. If the walls of the capillary are hydrophilic the 
liquid will move to fill it with a mean velocity (assuming 

two dimensional Poiseuille flow), v = —75^^, where ^ 
is the pressure gradient that sets up the flow. The driv- 
ing force for the flow is provided by the decrease in free 
energy as the fluid wets the walls or, equivalently, the 
Laplace pressure across the interface. Hence = — -rj, 
where 7 is the surface tension of the liquid interface, 
R — h/2cos9a is its radius of curvature and I is the 
length of liquid in the tube. Substituting ^ into the 
parabolic fiow profile and using v = dl/dt gives Wash- 



burn's law I — {'jhcosda/'ir])^^'^ {t + to)^^'^ , where to is 
an integration constant. Note that it is appropriate to 
use, not the static, but the advancing contact angle Oa, 
as this controls the curvature of the interface and hence 
the Laplace pressure. Washburn's law assumes that there 
is no resistance to motion from any fluid already in the 
capillary. The derivation also assumes incompressibility, 
and neglects inertial effects 01 , gravity f§\ and any de- 
viations from a Poiseuille flow proflle at the inlet or the 
interface 0,i,E,lii|. 

Our aim is to explore capillary filling in the presence 
of topological patterning on the surface of the capillary. 
Simple arguments, based on the Gibbs' criterion 5], will 
suffice to explain when capillary filling is not possible for 
one dimensional surface patterning (ridges). However, to 
see how the rate of filling is affected by the surface relief, 
and to consider two dimensional surface topologies, we 
will need to solve the hydrodynamic equations of motion 
of the fluids. Therefore we first describe the numerical 
model we shall use, a diffuse interface representation of 
a binary fluid, and check that it reproduces Washburn's 
law for simple capillary filling. 

We consider an immiscible binary fluid, comprising 
components A and B say, described, in equilibrium, by 
the Landau free energy 
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where p = pA + pB the total fluid density, and (f) = 
PA — PB is the concentration. We chose a — —(3 so that 
the two possible bulk concentrations are (f> — 1.0 and 
(f> = —1. 0. K is related to the surface tension by 7 = 
y^8Ka/9 and T is chosen to be 1 /3 to minimise the error 
terms. The surface contribution allows us to tune the 
equilibrium contact angle by varying a 



The fluid motion is described by the continuity, Navier- 
Stokes and convection-diffusion equations 

dtp + dppup = 0, (2) 
dtpua + dppuaup = -dpPap + dpS^p, (3) 

at<^(x, t) + do^^bUa, = Mdlpi^, t). (4) 
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FIG. 1: Lattice Boltzmann simulation results for capillary 
filling for h = 50, Oe = 60°, 7 = 0.0188, pA = ps = 1.0, 
rjA = 0.83, r/B = 0.03 and M = 0.05. (a) The length of the 
column of the filling A-component plotted against time. The 
circles are the simulation results and the solid line is a fit 
to the Washburn's law using the measured advancing contact 
angle and correcting for the small viscosity of the displaced 
B-component. (b) The measured advancing contact angle of 
the liquid-liquid interface: note that cos 9a varies linearly with 
capillary number [13I ]. 



where u is the fluid velocity, Pafs and Sap are the pressure 
and stress tensors respectively, fi is the chemical potential 
and M is a mobility. We solve these equations using a 
lattice Boltzmann algorithm, described in detail in [l2i |. 
Similar models have also been recently used in 0, Spto 
simulate capillary filling on smooth surfaces. 

Numerical results showing capillary filling of a smooth 
channel are presented in Fig. [TJa). The plot is for a 
channel of length L — 640, infinite width and height h — 
50. Reservoirs (480 x 200) of components A and B are 
attached at each end of the capillary. The two reservoirs 
are connected to ensure that they have the same pressure. 
The parameters of the model (given in the caption of Fig. 
[1]) are chosen so that the capillary and Reynolds numbers 
are of order 10~^. 

The solid line in Fig. [^a) is a fit to the Washburn's 
law using the measured values of the contact angle and 
correcting for the small viscosity of the displaced B- 
component. The fit is excellent, except very close to the 
beginning of the simulation where the expected devia- 
tions from Washburn's law due to inertial effects and de- 
viations from a Poiseuille flow profile are observed. Fig. 
[T](b) shows the measured cos 9a which decreases linearly 
as the capillary number decreases [isj . 

We now discuss capillary filling in a channel with topo- 
logical patterning on the surface. On stepped surfaces 
the driving force of the Laplace pressure must be suffi- 
cient to overcome any hysteresis due to the surface pat- 
terning. We shall assume that the contact angle of the 
A-component is < 90° i.e. this is the fiuid that may fill 
the capillary, and that there is no intrinsic hysteresis. 

We first consider two dimensions, and a channel where 



one of the walls is patterned by square posts, as shown 
in Fig. UK a). (The same results apply for a three di- 
mensional channel with ridges which are translationally 
invariant perpendicular to the flow.) According to the 
Gibbs' criterion, for an interface moving very slowly, the 
contact line will remain pinned on the edge of the posts 
until 9p, the angle the interface makes with the side of 
the posts, reaches the equilibrium contact angle 9(,. It 
follows immediately (see Fig. [SJa)) that there will be no 
capillary filling unless 

9p EE 90° ~9e< 9e, (5) 

that is, 9e < 45°. If the surface is patterned by triangular 
posts, with sides at an angle a to the edges of the channel 
as shown in Fig. Hfb) , the condition for capillarity is less 
stringent, namely 9^ < 90° — a/2. Hence the surface in 
Fig-Htc) will allow filling from left to right, but not from 
right to left if 45° < 9^ < 90° - a/2. 

Numerical results showing how the interface behaves 
as it approaches a post with finite velocity are presented 
in Fig. H^d) for channel height h = 50, post height 
hp — 10, post width Wp — 40 and equilibrium contact 
angles 9^ — 30° and 50°. Here I, the length of the liquid 
column, is measured at the centre of the capillary chan- 
nel. As the interface approaches the post, it is slowed 
down by dissipation caused by the post. However, once 
the interface is able to respond to the wetting potential 
of the post, it quicky wets the leading edge. The inter- 
face slows as it moves over the top the post ~ if the post 
was infinite in extent it would reach the speed appropri- 
ate to a channel of height h — hp. Once the trailing edge 
of the post is reached there is a substantial decrease in 
velocity as the interface becomes pinned. For 9e = 30° 
the interface is able to proceed, but for 9e = 50° it re- 
mains pinned, consistent with the condition in Eq. 
For the same value of 9e, the variation in velocity is more 
pronounced for higher posts. 

Fig. ^e) shows the dynamics as the interface moves 
along a capillary with groups of posts on both the top 
and bottom surfaces of the channel. This simulation 
demonstrates several features of the motion. First, as 
the interface reaches the three posts on the lower edge 
of the channel, it takes longer to depin from the first 
post it reaches. This is because the free end of the in- 
terface takes time to slide along the top of the capillary 
until it reaches an angle where it can pull the pinned end 
away from the trailing edge of the post. For the second 
and third posts the interface is already aligned diago- 
nally across the channel and does not need to move so 
far. Once the interface reaches the posts at the top of 
the channel it remains pinned on the first of these for 
some time until the interface again takes up a diagonal 
configuration, but now with the bottom edge leading. 

Such switching between two diagonally aligned inter- 
face configurations has been observed in experiments on 
capillary filling in channels with lateral root-mean-square 
roughness 1 /im, as shown in Fig. Hff). The chan- 
nels were realised by optical lithography employing plas- 
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FIG. 2: Pinning of the contact line during capillary filling when the top wall is patterned with (a) squares, (b) triangles, and 
(c) ratchets, (d) Interface velocity dl/dt when one of the walls is decorated with a single 40 x 10 post, (e) Interface velocity 
when the bottom wall is decorated with three posts and the top wall is decorated with two posts. The posts have height 10 
and width 5. (f) Experimental results on capillary filling in patterned channels. Local pinning temporarily holds the capillary 
front on one of the two microchannel sidewalls, while the interface advances along the opposite sidewall. 



tic photomasks patterned with a laser printer. Polymer 
microchannels were then obtained by subsequent rephca 
molding with polydimethilsiloxane. A 2 /il drop of water 
was released at the entrance of the microchannels and the 
penetrating interface was imaged by a 320x240 pixels, 30 
fps camera. Experiments have also shown the interface 
waiting longer at the beginning of a line of posts engi- 
neered to lie on one surface of a microchannel. 

We now consider two-dimensional surface patterning. 
Both walls of the capillary are patterned with 20 x 20 
square posts of height 10. The posts face each other 
across the channel, which is of height /i = 50. They are 
separated by spaces of length D in the direction perpen- 
dicular to the flow. This geometry is illustrated in Fig 
[3]Ja), which shows a side view of the channel. For Z? = 0, 
the case corresponding to ridges across the channel, the 
Gibbs' criterion predicts that the interface will always 
remain pinned between the trailing edges of a pair of op- 
posing posts. For D ^ 0, however, the motion can be 
driven by the portion of the interface that lies between 
the posts. 

Fig El^b) shows the crossover between capillary filling 
and contact line pinning as a function of D and O,,- For 
6'e < 30°, a small gap in between the posts is enough to 



allow filling. As the interface moves through the gap and 
protrudes beyond the posts, it is able to depin, initially 
from the side of the posts. 

For 6'e > 60°, however, filling is strongly suppressed. 
This is due to two complementary effects: the force driv- 
ing the interface though the gaps becomes weaker and 
the pinning on the sides of the posts becomes stronger. 
Note that the part of the interface that lies in between 
the posts can bow out further as D is increased and it is 
possible that the leading interface could touch the next 
row of posts, thus providing an alternative mechanism to 
depin the interface from the edges of the posts. 

We have also varied the relative height of the posts 
with respect to the channel width and find that the con- 
tact line pinning condition depends only very weakly on 
the channel height. This is because the interface be- 
haviour is primarily determined by depinning from the 
side edges of the posts. It is also worth noting that we 
neglect fluctuations in our simulations. 

So far, we have considered surfaces patterned with 
ridges oriented perpendicular to the flow and square 
posts, cases where pinning by the posts suppresses, or 
even prevents capillary filling. We now identify a situa- 
tion where patterning can be used to speed up capillary 
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FIG. 3: Two dimensional surface patterning: (a) the surface 
geometry; (b) regions of contact line pinning and capillary 
flow as a function of post separation and equilibrium contact 
angle. The channel height is kept constant at ft = 50. The 
dashed line is a guide to the eye. 




h = 50, M',, = 5. A,, = S,D = 10 
h = 50, = 8,h^, = 5,D = 7 
h = 50, w„ = 10. /i„ = 10D= 20 
h - 50, no patterning 
h - 40, no patterning 
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FIG. 4: The length of the column of ^-liquid plotted against 
time for different surface patterning. The channels' intrinsic 
contact angle is 9e = 30°. Other parameters as in Fig. [T] 



flow, ridges aligned along the flow direction. In Fig. IH we 
compare the lengths of the A column, plotted as a func- 
tion of time, for a channel of height /i = 50 with infinitely 



long ridges of width (wp), height {hp) and separation (D) 
(i) 5, 5, and 10, (ii) 8, 5, and 7, and (iii) 10, 10, and 20 ly- 
ing along the channel. In all these cases the channels have 
the same surface roughness r — \ + 2hp/{wp + D) — 5/3. 
For comparison, we also show smooth channel simula- 
tions for h = 40 and h = 50. It is clear from Fig. [H that 
it is possible to increase the rate of capillary filling by 
patterning the channels (cases (i) and (ii)). This occurs 
because there is no contact angle hysteresis in this geom- 
etry and because the total surface area per unit length 
wetted by the ^-liquid is increased over the flat channel 
value by the roughness factor r. However, from Fig. [Jit 
is apparent that the rate of filling is not just controlled by 
r. This is because surface roughness distorts the liquid 
fiow from a parabolic profile and hence increases the dis- 
sipation. As shown in Fig. [51 the dissipation increases as 
the distance between two neighbouring ridges is reduced 
(case (ii)) and as the size of the ridges relative to the 
channel height is increased (case (iii)). The filling speed 
is a competition between the increased energy gain of 
wetting a grooved channel and the additional dissipation 
caused by the distortion of the fiow. 

Our results suggest strategies that could be used to 
design surface patterning on microfluidics channels. To 
speed up capillary filling, one should introduce roughness 
which is elongated in the direction of the flow. However, 
additional surface area will only speed up filling if pin- 
ning is avoided by, for example, avoiding sharp edges and 
making the post inclinations as gentle as possible. Fur- 
thermore, any surface patterning should neither be too 
dense nor too large to minimise dissipation. Random 
arrangement of the posts may also be used to avoid co- 
operative interface pinning effect. Finally, the intrinsic 
contact angle of the materials 6^ < 60° for capillary fill- 
ing down a rough channel. These results are consistent 
with recent findings by Kohonen [T^ . where he found 
that the walls of the tracheids in the dry-habitat species 
are typically rougher than those in wet-habitat species, 
and with 9^ ~ 40°. Moreover, the set of design principles 
we describe above may help explain why certain types of 
wall sculpturing, such as helical thickening found on 
tree capillaries, are more common than others. 
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